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ABSTRACT 

A modeling approach for linear isentropic 
flow systems based on the quasi-one dimensional 
Euler equations of non-viscous, compressible flow 
is presented. Such systems are representative of 
certain high speed propulsion systems. Accurate 
models useful in control system studies are 
developed. A supersonic inlet is considered and the 
resulting set of partial differential equations with 
boundary conditions is solved for a linear transfer 
matrix using Laplace transforms. 


INTRODUCTION 

High speed propulsion systems, such as 
supersonic inlets, can be described by the quasi- 
one dimensional set of Euler equations which 
represent non-viscous compressible flow. Quasi- 
one dimensional flow occurs in a channel with a 
small rate of change of cross section (see Figure 1 
beiow) or between nearly parallel streamlines in 
which case the velocity components normal to the 
mean flow direction are small compared to the 
components parallel to the flow direction. 



Figure 1. A Stream Tube. 


The resulting set of equations is nonlinear 
and includes equations for conservation of mass, 
conservation of momentum, and conservation of 
energy. Usually, any solution to this set of 
equations requires computational fluid dynamics 
(CFD). However, if it is assumed that any 
perturbation in the flow field is small, a linearized 
set of equations results. Furthermore, if it is also 
assumed that any energy wave is insignificant, then 
the conservation of energy equation can be 
eliminated from the set. This situation is called 
isentropic flow ; as the entropy remains constant in 
the flow field. 


The following linear isentropic flow equations now 
result: 

at ax p a x 

i£. + p§ii + u i>£ = f p 

at ax ax 

where u is the flow velocity, p is the density of the 
medium, the fs are spatially distributed inputs, c 
is the speed of sound, and the barred variables 
represent the steady state values around which 
these equations were linearized, 

BVP PROBLEM FOR A SUPERSONIC INLET 


To apply the linear isentropic flow 
equations to a real system, boundary conditions are 
necessary. For a specific solution, a supersonic 
inlet is considered. Here the downstream boundary 
is simply reflective. In order to represent 
compressor or combustor signals, a perturbation 
input f ds is also included in the downstream 
boundary. On the other hand, the upstream 
boundary is considered to be the normal shock 
usually located near the throat of a supersonic 
inlet. Furthermore, to represent perturbations in 
the free stream atmosphere, an upstream input f U5 
is included. Symbolically, these conditions are 

Pds=R Uds+fds» Uu$=Zus(s) p us + fus, 
where R is a constant representing the reflective 
boundary and Zus=Zus(s) is a lumped dynamic 
representing the shock behavior (see Figure 2). 

Using zero initial conditions and Laplace- 
transforming the given system yields: 

S u(s,x) + uiHii^<£Mi^) = f u 

9x p ax 

sp(s.x)+p* u ( s Au?P( s - x )=f,, 
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If f u =f p =0, then using matrix notation the above 
system becomes, 

A^S^2 + By( S ,x)=o 

0 X 

where, 


A = 


L P 


P 
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Since, matrix A is nonsingular, the system can be 
rewritten as, 

9 4 s - Vbks, X) = 0 

3 x 

The boundary conditions for this system, using the 
notation of reference 1, can be put in the following 
form, 

M(s) Q(s, 0) + N(s) Q(s, L) = Z(s) 
or, 

M(s) Qu. + N(s) Q ds - Z(s) 

where, 



Note that since Zfs^O, the termination is said to 
be "nonhomogeneous." The solution of the system is 
u(s, x) = e' A ’ 8 * u(s, 0)3 H(s, x) y(s, 0)3 H(s, x) yu. 


1 -Z u , 
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Thus, the transfer function matrix K(s,x) is given 

by. 

K(s,x)=H(s,x)[M(s) + N(s)H(s,L)] ’’eKo Ki Kz 


Ka K 4 J 


where, the quantities K 0 , K tl K 2 , K 3 , and K 4 are 
given next. 
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where, the open loop or "free space" transfer 
function matrix H(s,x) is given by 
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The existence of hyperbolic functions in H(s,x) is 
an indication that the resulting model is a 
complicated transfer function having an infinite 
number of poles and zeros, characteristic of 
distributed parameter systems. The transfer 
function matrix from Z(s) to Q(s, x) can be shown 
to be given by 

K(s, x) = H(s, x) [M(s) + N(s) H(s, L)] 1 
Thus, calculation of K(s,x) requires first 
computing H(SsL). From H(s,x), given above, it 
follows that. 


H(s, L)=ex 




Furthermore, substituting M(s), N(s), and H(s,L) 
into K(s,x) one obtains the following expressions. 
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inally, the bounded system is as follows: 
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K(s,x) is referred to as the transfer function 
matrix for the given boundary conditions and is 
given by, 


K(s, x)« 


Ku Ki 2 


L K 2 , K22 



It is interesting to note that every term in K(s,x) 
is a ratio of entire functions. This is true because 
the quantity M(s)+N(s)H(s,L) as well as its 
adjoint is a sum of entire functions. 


NOISE 



Figure 2. A Supersonic Inlet System. 


APPLICATION 

From the previous section, it can be seen 
that Q(s,x) is the output vector of the closed loop 
system when its input is the vector 

Z(s)=^fus fds] 


f ug and f dg will be assumed to be two noise 
corrupted control inputs to the supersonic inlet 
system. This situation is illustrated in Figure 2. 

it should be noted that the system transfer 
function K(s,x) can be more complicated than it 
appears since Zu$ is a lumped dynamic system 
representing the shock behavior. In general, this 
can be a ratio of two polynomials having m zeros 
and n poles. The corresponding mathematical model 
for Zus would then be as follows: 


Zus(s) = 
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where, it is usually assumed that ao^O, and bo*0. 
However, it has been shown experimentally (see 
references 2 and 3) that Zus(s) can be modeled as a 
first order system having a fast lag. This suggests 
that Zus(s) takes the following form: 


Using the data shown in Table 1, which can be 
obtained from reference 4, one can compute the 
frequency response of the closed loop system at 
x=0, which will give an approximate response for 
the shock. 


Table 1. Supersonic Inlet Data 


u 

574 

p 

0.00136 

c 

113 5 i 

L 

2.57 


The shock transfer function is assumed to be 

Zus(S) - IQ 6 

s + 3 0 0 

and the reflection parameter, R, is chosen to be 
R=6-10' 7 . At this point, the values for R and Zo 
are roughly chosen so that the magnitude of the 
reflections* from the boundaries are approximately 
unity. Further research to develop a more 
accurate technique for determining R and Z 0 as 
well as better shock models is currently being 
conducted. 

Using the above data the response shown in 
Figure 3 is obtained. Note that only the K 12 (s,0) 
term is shown. For comparison purposes the 
frequency response of this system computed in 
reference 4 is superimposed on the frequency 
response of Ki2(s,0). 

Furthermore Figure 4 illustrates a three- 
dimensional plot of log|Ki2(s,0)| versus s. The 
peaks in the plot correspond to the pole locations of 
Ki 2 (s,0). 
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Figure 3 . Magnitude Response for IK^s.O}). 



CONCLUSIONS 

The resulting model is very useful to the 
control engineer. It can provide quick fast insight 
about this very complicated system. The solution 
obtained here is in a compact form and one can get 
the answer quickly. In fact, a MATLAB program 
has been developed to compute frequency responses 
for various conditions and shock models, which are 
in turn useful for control system design using root 
locus or frequency domain techniques. 

The analytical solution obtained here is an 
exact solution of the boundary value problem as 
opposed to the approximate solutions presented by 
Cole and Willoh in references 2 and 3. The solution 
is good as long as the conditions for a small 
perturbation model are met. However, to obtain a 
better model, the energy equation should be added to 
the set. 

Finally, it is noted that the free space 
system function matrix given in this paper is very 
similar to that of an electrical transmission line 
system. In fact, both the inlet and the transmission 
line systems are in a sense identical for u = 0. 
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Figure 4. loglK^s.O)! versus s. 

To better understand the location of the poles, a 
contour plot is also constructed in Figure 5. Closed 
curves represent points in the s plane that yield 
the same magnitude in the frequency response. 
Note that all the poles seem to lie in a strip 
approximately between s=-150 and s=-350 and 
they are symmetrical about the real axis. 
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